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Abstract 

We compute the tail asymptotics of the product of a beta random variable and a generalized 
gamma random variable which are independent and have general parameters. A special case of 
these asymptotics were proved and used in a recent work of Bubeck, Mossel, and Racz in order 
to determine the tail asymptotics of the maximum degree of the preferential attachment tree. 
The proof presented here is simpler and highlights why these asymptotics hold. 


1 Introduction 

There has been a lot of recent interest in various urn schemes due to their appearance in many graph 
growth models (see, e.g., [isi m o da n da 0 El IS]). The limiting distributions arising in these 
urn schemes are often related to the beta and gamma distributions. Consequently, the computation 
of various statistics in random graph models often boils down to using algebraic properties of these 
distributions, commonly referred to as the beta-gamma algebra [6]- 

The purpose of this note is to simplify and demystify a recent computation done in [5] involving 
beta and generalized gamma random variables. Bubeck, Mossel, and Racz [5] were interested in 
the influence of the seed graph in the preferential attachment model, which led them to study 
the tail asymptotics of the maximum degree of the preferential attachment tree. This, in turn, 
essentially reduces to computing the asymptotics oiF{BZ > t) as t ^ oo, where B ~ Beta (o, 6) 
and Z ~ GGa (a + 6 + 1,2) are independent random variables and a and b are positive integers; 
here Beta (a, b) denotes the beta distribution with positive parameters a and b (with density 
B(a where B{a,b) = is the beta function), and GGa(c,p) de¬ 
notes the generalized gamma distribution with density for c, p > 0. We refer 

to |5j for dots-ils on tlioso connoctionsj soo nlso [iniiiii. 

The computation in [5| involves a few pages of alternating sums cancelling each other out in 
just the right way. Here, in contrast, we provide a short and simple proof of these asymptotics. 
The core calculation is only a few lines long, involving approximations at three points which are 
natural and which can be justified in a relatively straightforward manner. Moreover, the argument 
works for all positive values of the parameters a, b, c, and p. Throughout the paper we use standard 
asymptotic notation; for instance, f {t) ^ g {t) as t — >■ oo if limt^oo f (t) /g (t) = 1. 

Claim 1. Let a, b, c, and p be positive, let B ~ Beta (a, b), and let Z ~ GGa (c,p), with B and Z 
independent. Then we have 

P {BZ > t) ~ ast^oo. 

T (c/p) r (a) 
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There has been lots of work on understanding the distribution and tail asymptotics of products 
of random variables; see, e.g., [I6] for a paper from nearly half a century ago, and [9] and references 
therein for recent developments. In particular. Claim [1] is a special case of O Theorem 4.1], where 
the authors prove a general result for any product BZ where B ~ Beta (a, b) and Z has a law which 
is in the maximum domain of attraction of the Gumbel distribution. Due to the generality of their 
result their proof is fairly involved. We thus believe that the simple proof we present here is useful 
in highlighting why these asymptotics hold. 

The product BZ studied in Claim [T] has many nice properties, for instance it has moments 
of Gamma type m- When p = 1, Z is a gamma random variable and BZ has a so-called G 
distribution, with its moments described by Meijer’s G-function [7]. The case of p = 2 appears in 
many settings, including the preferential attachment model as mentioned above, and see also mm 
for connections to critical random graphs, random walks, and various random trees, including 
Aldous’s Brownian continuum random tree (CRT). In a very interesting recent work, Pekoz, Rollin, 
and Ross m showed that generalized gamma random variables with p being an integer greater 
than 2 arise as limits in time inhomogeneous Polya-type urn schemes. It would be interesting to 
find connections to urn schemes and random graph models for general values of p. 


2 The core calculation 


In this section we prove Claim [T] modulo some approximations whose validity is justified later in 
Section [3l Put W = Z^, z = and dz = The density of W is thus 

In other words, W has a gamma distribution: W ~ Gamma (c/p, 1). We have 

P {BZ >t)=¥ {BPZP >tP)=¥ {BPW > tP ), 


and so the claim is equivalent to 

P {BPfV > tc) ~ y ^ p~^w‘^^P~^~^e~'^ as w ^ oo. 

T {c/p) T {a) 

We can write 

P {BPW > w) = F{W > w)F {BPW >w\W>w). 

The first factor has well-known asymptotics: P {W > tc) ~ fw {w) as rc —>■ oo [H formula 6.5.32]. 
Also, it is well known that for any random variable W which has a Gamma (r, 1) distribution, we 
have 

{W — w\W > w) ^ £ as w ^ oo. 


where T is a standard exponential random variable. This convergence is rather strong, e.g., con¬ 
vergence of densities. So 


P {BPW > w\W > w) = F {BP {w + W - w) > w\W > w) 


F{BP{w + £)>w)=f{b> 


1 


F{ B >1- 


pw 


= P 1 - R < 


1 + £/w 

£ 


i/p^ 


pw 


poo 1 1 / 1 \ ^ 

/ rTty — r(6) 

Jo B{a,b) B{a,b) \pwJ 


( 1 ) 

( 2 ) 

(3) 

(4) 
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and the conclusion follows. The only thing that remains is to rigorously justify the three points 
where asymptotic equivalence was used in the line of reasoning above; see Section [3] for details. 


3 Justifying the approximations 


Here we justify why the expressions in ([T]), ([2]), ([3]), and ([4]) are all asymptotically equivalent as 
w —)• oo. 

Asymptotic equivalence of ([3|) and Q. By definition we have 


P ( 1 - H < —) = [ ] 

pwj Jo 


1-B < — ] e-^dz = 
pw 


I ' 


— ] e-^dz + e-P^. 
pw ) 


Since u^~^e~P‘^^du = 0 (w~^) as w ^ oo, we can neglect all terms that are o [w~^) as w ^ oo, 
and so we only have to deal with the integral term on the right hand side of the display above. 
Using that 1 — B ~ Beta (6, a) and also the change of variables z = pwu, we have 


ppw 

I ' 


l-B < — ] e-^dz = 
pw 


pw p 1 

/ P'^ J. 


ppw P 

Jo Jo 


^ (1 — x)°' ^ dxe ^dz 


a—1 


pw 


B (a, b) 

1 pU 



x’’-^ (1 - dxe-P^^du. 


B {a, b) Jo Jo 

The integral from 1/2 to 1 is negligible, since 


pw 


t — 

Ji/2 B (a. 


x^ ^ (1 — xY 


b) 


Uxe-P^^du< C {pw)e-P^'^du<e-P^/‘^. 
Jll2 


For the integral from 0 to 1/2 we first drop the factor (1 — x)“ ^ and we justify the validity of this 
later. Integration by parts then tells us that 


"1/2 


^{pw)e-P^^du =+ / vJ’-^e-P^^du 
0 2 b Jq 


poo 

/ 

Jo 


e-P^^du - -^e-P^/'^ - 
2’^b 


poo 

/ 

4i/2 


e-P^^du. 


Using that (1/2 + ^ < 2 + (2n)^ ^ for all b > 0 and n > 0, we have that the integral from 1/2 

to oo in the display above is negligible: 


POO 

/ u'^-^e-P'^'^du = 

/l/2 


, 6-1 


POO 

/ (1/2 + ^)' 

— piu /2 i: {2 + ( 2 n)^-^) e-P'^^dv < poly {w) , 


-pwl2 


where poly {w) is some polynomial in w. We have thus shown that 

'* 1/2 pu poo 


/ r x’^-^dxe-P^^du^ r / u^-^e-P^'^du 

B{a,b)Jo Jo Jo B{a,b) 


as w —> 00 . 
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Finally, to justify dropping the (1 — x)“ ^ factor, note that 1 — (1 — x)“ < max {a, 2} x for all 

a > 0 and x G (0,1/2), and so 


pw 


^1/2 pu 

/ / x^-^ (1 - xf-^ dxe-P^^du - 

Jo Jo 


B (a, b) 


pw 


" 1/2 nu 


10 JO 
r-1/2 pu 


pu 

/ x'^-^ 

Jo 


dxe-P^^du 


< max 


B (a, b) 

p \_ j p Uj 

{a, 2} / / x^dxe-P^^du = O . 

B [a, b) Jo Jo V y 


Asymptotic equivalence of ([2|) and ([3]). We need to show that, as w ^ oo. 


B > 1- — ] -F \ B > 
pw 


1 + £/w 


l/pN 


= o \w 


..-b 


(5) 


Using that 1 — x/p < (1 + x) < 1 — x/p + x > 0, we have that the left hand side 

of ([5]) is at most 


B G 


1-, 1-h 


£ p+l £'^ 


pw ’ pw 2p^ w"^ 


,£ < ) + 2 f(£> -^w 

p+i j V p+^ 


We know that P > -pppwj = e p+i“^, so we only have to deal with the first term in the sum 
above. Using the change of variables z = pwu we have 


P UB G 


£ ^ £ p + 


1-, 1-h 


pw ’ pw 2p^ w"^ 


,£< 


p 


P + i 


-w 


P p — 

j P +1 / Pt 

Jo J ^ 


1 


x^ ^ (1 — x)“ ^ dxe ^dz 


a—1 


-E±^zl B {a, b) 

pw 2 p‘^ ^ 

1 


pw 


r x^’-^ (1 - x)“-^ dxe-P"'“dn 

Jo Ju-E^u‘2 


< 


B (a, b) 

{l,(p/p+l)“-^} 


max 


B (a, b) 


{pw) f 

Jo J u 


x^-^dxe-P^'^du. 


0 


Using simple estimates we have = i for all 


e ( Or ^ j and so 


{pw) 


. 1 
p+1 


pu 

Ju-E^i 


x^-^dxe-P^'^du < 


10 


< 


(p + 1) (6 + 1) 

2h 


1 

p+1 


{pw) / vJ'+^e-P^^du 


(p+1) (6 + 1) 

2b 


{pw) u’^+^e-P-^du = 


which concludes the proof of (|5]l. 

Asymptotic equivalence of ([T]) and ([2]). We justify this via a direct calculation, though 
there might be a more elegant way to do this. To abbreviate notation, let r := c/p. By definition 
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we have 


¥ {BP {w + W - w) > w\W > w) 


¥{BP{w + y-w)> w) Pdy 

,oc ; - r ¥{BP{w + z)>w)(l + -) e-^dz. 

x^-^e-^dx Jo wJ 


For r = 1 this is exactly equal to P {B^ {w + £) > w). Since P {W > tc) ~ fw (w) as w ^ oo, the 
fraction in front of the integral in the display above goes to 1 as tc ^ oo, so what remains to show 
is that 


P (BP {w + z) > w) 


1 + 


w 


r—1 


- 1 


e ^dz = o (P {BP {w + £) > w)) 


as w 


oo. 


We partition the integral into two parts: from 0 to ^/w, and from ^/w to oo. For the first term 
note that for r > 1 and z G (0, ^/w) we have (1 + z/wY~^ — 1 < — 1, while for r G (0,1) and 

z G (0, y/w) we have 1 — (1 + zlwY~^ < 1 — (1 + l/^/w)~^ < l|^/w, and so 


rypw 


/ P {BP {w + z) > w)\ 

(l + -) -1 

/o 

V w/ 


e ^dz 


< max 


l| J F {BP {w + z) > w) e ^dz 

< max < -^=,e^ — 11 P {BP {w + £) > w) = o {P {BP {w + £) > w)) . 
I Vrc J 


For the second term we can bound the factor P {BP {w + z) > w) in the integral by 1. For r G (0,1) 
we have 1 — (1 + z/wY~^ < 1 — (1 + z/w)~^ < z/w and so we get the upper bound of 


-j /•OO -j /•OO 

- / ze-^dz = -e-^ / {^/^ + u) e-^du = 

w w Jo '■ 


tc + 1 _ 


w 


For r > 1 we use the bound {1 + z/wY ^ — to obtain the upper bound of 

/•OO 1 

/ Y^/^-^^^dz = - < 2e-v^/^ 

where we assumed that w > 2r. 
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